Fine structure of charge-exchange spin-dipole excitations in ie O 
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The charge-exchange spin-dipole (SD) excitations for both (p, n) and (n,p) channels in ls O are 
investigated in the fully self-consistent random phase approximation based on the covariant density 
functional theory. The fine structure of SD excitations in the most up-to-date 16 0(p, n) 16 F experi- 
ment is excellently reproduced without any readjustment in the functional. The SD excitations are 
characterized by the delicate balance between the a- and o;-meson fields via the exchange terms. 
The fine structure of SD excitations for the le O(n,p) 16 N channel is predicted for future experiments. 
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The nuclear charge-exchange excitations [l[ correspond 
to the transitions from the ground state of the nu- 
cleus (N, Z) to the final states in the neighboring nuclei 
(N =p 1, Z ± 1) in the isospin lowering T_ and raising T + 
channels, respectively. These excitations can take place 
spontaneously such as the well-known (3 decays or be 
induced by external fields such as the charge-exchange 
(p, n) or (n,p) reactions. They are categorized accord- 
ing to the orbital angular momentum transfer as allowed 
transitions with L = and first- and second-forbidden 
transitions with L = 1 and L = 2, etc. Meanwhile, they 
are also classified by the spin's degree of freedom as the 
non-spin-flip modes with 5 = and the spin-flip modes 
with S = 1. 

Among all the nuclear charge- exchange excitation 
modes, the spin-dipole (SD) excitations with S = 1 and 
L = 1 have attracted more and more attentions due to 
their connection with the neutron-skin thickness I2J, the 
cross sections of neutrino- nucleus scattering 0, 0], the 
double beta decay rates Q , and so on. Different from the 
famous Gamow- Teller (GT) excitations having a single 
spin-parity J 71 = 1 + component, the SD excitations are 
composed of three collective components with spin-parity 
J n = 0~,1 _ , and 2~. It is relatively straight forward 
to distinguish the orbital angular momentum transfer L 
by the angular distributions of double differential cross 
sections, but it is not trivial to resolve different J n com- 
ponents in SD excitations [l|. However, the resolution 
of these three J T components is crucial to understand 
the multipole-dependcnt effects on the neutrino-nucleus 
scattering Q and neutrinoless double beta (OV/3/3) de- 
cays @, HI , and the strengths of nucleon-nucleon effective 
tensor interactions for understanding the evolution 

of the single-particle energies in exotic nuclei 11 , 12| ■ In 
particular, the J T = 0~ states can also serve as door- 
ways for parity mixing in compound nuclear states fl3j |. 
Therefore, the investigation of the fine structure for SD 
excitations including all components has become one 
of the central issues for both experimental and theoretical 



nuclear physics, particle physics, and astrophysics. 

Charge-exchange excitations in 16 are of particular 
interest in both nuclear physics and astrophysics. For 
instance, 16 is the key nucleus in the water Cherenkov 
detector for (anti-) neutrinos providing evident signals of 
supernova neutrino bursts and neutrino oscillations 0- 
1 1 61 ] - In one of the most recent experiments, the fine struc- 
ture of GT and SD excitations in 16 has been identi- 
fied by using the 16 0(p, n) 16 F reaction with polarized 
proton beam 171 ]. In this experiment, the known SD 
< 8 MeV have been clearly identified, 
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where the excitation energies E x are measured from the 
ground state of the daughter nucleus 16 F. As shown with 
arrows in Fig. [TJ the peak at E x ss MeV is composed 
of the triplets of J 71 ' = _ , 1~, 2~ states, while the main 
SD resonance at E x w 7.5 MeV and the "shoulder" at 
E x = 5.86 MeV are found to be J n = 2~ states. It is 
also identified from this experiment that the broad SD 
resonances at E x 9.5 and 12 MeV are formed by the 
mixture of J x = 1~ and 2~ states, where the former 
resonance is dominated by the J 11 — 2~ component and 
the latter one is dominated by the J v = 1~ component. 
The experimental data in such details provide a rigorous 
calibration for theories, in particular the microscopic the- 
ories aiming at describing both ground-states and excited 
states all over the periodic table with a high predictive 
power. 

Theoretically, the nuclear charge-exchange excitations 
are mainly investigated by the shell model and the ran- 
dom phase approximation (RPA) built on energy density 
functionals. Limited by the computational facilities avail- 
able, the RPA approach is the only microscopic method 
that can be implemented for the whole nuclear chart. 
The same energy density functional should be used for de- 
scribing both the nuclear ground state and excited states 
for the model self-consistency [ljj [2(j, in order to re- 
store the symmetries of the system Hamiltonian broken 
by the mean-field approximation, to separate the spuri- 
ous states from the physical states, as well as to maintain 
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the predictive power for nuclei far away from the sta- 
bility line. Recently, such full self-consistency has been 
achieved in the framework of covariant density functional 
theory (CDFT) [2l|, i.e., the self-consistent RPA built on 
the relativistic Hartree-Fock (RHF) theory, denoting as 
RHF+RPA below. Excellent agreement with the GT res- 
onances data in 48 Ca, 90 Zr, and 208 Pb has been obtained 
without an y re adjustment of the covariant energy density 
functional (21| . 

In this paper, the self-consistent RHF+RPA approach 
will be used to investigate the fine structure of the SD 
excitations for both (p, n) and (n,p) channels in 16 with 
the most up-to-date data. The main focus will be the un- 
derstanding of the characteristics of SD excitations and 
the calibration of theoretical models. 

The basic ansatz of RHF theory is an effective La- 
grangian density £, in which nucleons are described as 
Dirac spinors that interact with each other by exchang- 
ing the a, u>, 
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The 



p, it mesons, and photons 
system energy functional E is then obtained as the ex- 
pectation of the effective Hamiltonian sandwiched by the 
trial ground-state wave function within Hartree-Fock and 
no-sea approximations. This theory is also called Dirac- 
Fock method in atomic physics [2; 



25|]. For the fully 



self-consistent RPA established beyond, the particle-hole 
(ph) residual interactions are strictly derived by taking 
the second derivative of the same energy functional E 

as fum, 



> 
s 



> 
s 

a. 




5 10 
E x (MeV) 



15 



20 




K(l,2) = -[ff«Oto]i[fe7o]2Ar(l,2), 

K,(l,2) = [^7071 1^707^2^(1, 2), 

K„(l,2) = [^707^1^707^2^(1,2), 

14(1,2) = -[I^T lGlbl k d k ]i 



T70757 i ^]2-D7r(l,2), 



(la) 
(lb) 
(1c) 



(Id) 



where -0,(1,2) denotes the Yukawa propagator. The pi- 
onic zero-range counter-term which cancels the contact 
interaction of the pseudovector ir-N coupling is 

14,5(1,2) = g'[— T7o7s7]i-[— T7o757]25(ri-r2), (2) 

where g' = 1/3. 

In the present RHF+RPA calculations, the effective 
interaction is taken as PKOl J23|, which is determined 
by a set of selected nuclear ground-state properties, and 
has no free parameters for nuclear excitations. The radial 
Dirac equations are solved in coordinate space within a 
spherical box with radius R = 25 fm and mesh size dr = 
0.1 fm 2jj. The single-particle energy truncations for 
RPA matrices are taken as [-M, M + 100 MeV] with M 
the nucleon mass. 

In the upper panel of Fig. [TJ the strength distributions 
of the SD excitations in the T_ channel of 16 O calculated 



FIG. 1: (Color online) Strength distributions of the SD exci- 
tations in 10 O for the T_ channel. The J n =0~, l"j_and 2" 
contributions calculated by RHF+RPA with PKOl (up- 
per panel) and RH+RPA with DD-ME2 [IH] (lower panel) 
are shown as the dash-dotted, dotted, and dashed lines, re- 
spectively, while their sums are shown as the solid lines. The 
energy of the lowest RPA state is taken as reference and a 
Lorentzian smearing parameter T = 1 MeV is used. The ex- 
perimental data [17I. [la] are shown with arrows, whose widths 
illustrate the widths of the corresponding resonances. 



by the RHF+RPA approach with PKOl j2J are shown 
by taking the lowest RPA state as reference. The spin- 
parity J 71 " = 0+ 1+ and 2~ contributions are shown as 
the dash-dotted, dotted, dashed lines, respectively, while 
their sum is shown as the solid line. For comparison, 
the experimental low-lying SD excitations [18l ] and reso- 
nances 17| are denoted with arrows, while their widths 
illustrate the widths of the SD resonances at E x w 7.5, 
9.5, and 12 MeV, respectively One should be careful 



that, different from the GT case, the strength distribu- 
tions are not proportional to the experimental cross sec- 
tions at zero-degrees for the SD excitations. 

In general, the profiles of the J n = 0~, 1~, and 2~ 
excitations are well reproduced by the calculations. Fo- 
cusing on the details, firstly, the 07, If, and 27 triplets 
are found at E x m MeV. Secondly, the main giant reso- 
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nance at E x ps 7.5 MeV as well as its "shoulder" structure 
at E x ps 6 MeV (which could not be described by the shell 
model calculations (T3|) generated by the J T = 2 _ com- 
ponent are excellently reproduced. The "shoulder" struc- 
ture is formed by the coherent excitations of (vp^hn d^/2) 

and iypy^d^j^)- In contrast, it is suppressed by the in- 
terference between iyp~^J 2 T:d^/2) and (yp~^^KS\j2) in the 
shell model calculations. Thirdly, the broad resonances 
at E x ps 9.5 and 11 ~ 13 MeV are understood as the 
mixture of the J 71 " = 1~ and 2~ excitations. Accord- 
ing to the transition strengths, the former resonance is 
dominated by the J 71 " = 2~ component, whereas the lat- 
ter one is dominated by the = 1~ component. To 
conclude, the fine structure of the SD excitations can be 
described robustly by the self-consistent RPA based on 
a covariant density functional without any readjustment. 
For the J v = 0~ resonances beyond E x = 10 MeV, where 
no clear bumps have been observed experimentally [l7| . 
the present theory predicts such resonances being frag- 
mented at the range of E x = 12 ~ 18 MeV with the peak 
at E x ps 14.5 MeV. The peak e nerg y is consistent with 
the shell model prediction in Ref. |l7J . Such a fragmented 
feature could be considered as one of the reasons why the 
J n = 0~ resonances are so difficult to be observed. 

To illustrate the effects of the exchange terms, the cal- 
culated SD strength distributions by the conventional 
RH+RPA approach [H with DD-ME2 [H] are shown 
in the lower panel of Fig. Q] where g' in Eq. ([2]) is read- 
justed to 0.52 according to Ref. [30[. It can be seen that 
the general profile of E x < 8 MeV calculated without 
exchange terms is similar to that with exchange terms. 
However, substantial discrepancies can be seen for the 
SD resonances beyond E x = 8 MeV. The mixture of 
the = 1~ and 2~ excitations at E x ps 9.5 MeV is 
missing, and the J 11 = 1~ resonances are too high in en- 
ergy by comparing with the data. It is also worthwhile 
to note that the J* = 0~ resonances are centralized at 
E x = 10 ~ 12 MeV. 

In order to understand the differences between the pre- 
dictions by RHF+RPA and RH+RPA, the unperturbed 
and collective SD excitations are shown in Fig. [2j The 
unperturbed excitations for all J 71 components are found 
to be quite similar by these two approaches. This indi- 
cates that there is no substantial difference for the cal- 
culated single-particle spectra with or without exchange 
terms. However, when the collectivity is switched on, it is 
clear that the ph residual interactions in these two frame- 
works play very different roles. In particular, the most 
profound difference is found in the J™ — 0~ component, 
i.e., the residual interactions are repulsive in RHF+RPA 
and slightly attractive in RH+RPA. 

To further understand these discrepancies and evaluate 
the significance of the model self-consistency, by taking 
the diagonal matrix elements of the main ph configura- 
tions for the J 71 = Q~ , 1~, and 2~ excitations as cxam- 
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FIG. 2: (Color online) The unperturbed (Hartree-Fock for 
PKOl and Hartree for DD-ME2) and collective (RPA) exci- 
tations for the J n — 0~, 1 _ , and 2~ components. 
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FIG. 3: (Color online) Diagonal matrix elements of the 
particle-hole residual interactions for the J 71 — 0~, 1~, and 
2~ excitations. The total strengths are decomposed into the 
contributions from the a and ui mesons (a+uj), the p me- 
son, the pseudovector n-N coupling (tt PV) and its zero-range 
counter-term (n ZR). The configurations are displayed as the 
format of neutron-hole-proton-particle. 
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pies, the total strengths of the ph residual interactions 
are presented in Fig. [3] together with the corresponding 
contributions from the a and lu mesons, the p meson, the 
pscudo vector ?r-N coupling and its zero-range counter- 
term. 

Within RHF+RPA, the total strengths of the ph resid- 
ual interactions are essentially determined by the deli- 
cate balance between the a and ui mesons via the ex- 
change terms. The p-meson seems to be important in 
the natural-parity channel, while the pseudovector Tt-N 
coupling and its zero-range counter-term play minor roles 
due to the strong suppression of the coupling strength 
of PKOl in the framework of density-dependent RHF [23] 
in nuclear medium. As the tr and u) mesons are well cal- 
ibrated by the central and spin-orbit potentials for the 
nuclear ground-state properties, the proper description 
of the SD excitations provides a stringent and critical 
test of the theoretical model. It is noted that, different 
from the Skyrme HF+RPA calculations , here the ex- 
plicit tensor interactions are not necessary to reproduce 
the data. 

Within RH+RPA, the total strengths of the ph resid- 
ual interactions are essentially determined by the p me- 
son for the natural-parity channel and by the pseudovec- 
tor ir-N coupling and its zero-range counter-term for the 
unnatural-parity channels. As the pseudovector ir-N cou- 
pling and its zero-range counter-term are absent in the 
description of the nuclear ground-state properties, the 
adjustment of <?' with the GT (J^ = 1 + ) excitation en- 
ergy introduces an extra parameter and therefore influ- 
ences its self-consistency. 

Encouraged by the success in reproducing the fine 
structure of charge-exchange SD excitations in the T_ 
channel of 16 O, it is worthwhile to investigate other 
phenomena interesting for future experiments by using 
the microscopic and fully self-consistent RHF+PRA ap- 
proach. 

The strength distributions of the SD excitations in the 
T + channel of 16 O calculated by the RHF+RPA with 
PKOl are shown in Fig. 0] by taking the lowest RPA 
state as reference. The 07, 17; an d 27 triplets are found 
at E x = ~ 2 MeV. The main resonance and a shoulder- 
like structure dominated by the J T = 2~ component ap- 
pear at E x ss 8.3 MeV and E x w 6.5 MeV, respectively. 
The broad giant resonance at E x = H ~ 15 MeV is su- 
perposed of the J T = 1~ and 2~ excitations with the 
J w = 1 _ component dominant in transition strengths. 
At E x = 13 ~ 17 MeV, the J" = CT resonances are 
predicted to be fragmented. These calculations are well 
supporte d by the available ex per imental low-lying SD ex- 
citations [18| and resonances |31[ denoted by arrows with 
their corresponding widths. The predicated fine struc- 
ture might be verified by future experiments with polar- 
ized beams. 

In summary, the charge-exchange SD excitations in 
16 O have been investigated with the fully self-consistent 
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FIG. 4: (Color online) Same as the upper panel of Fig.[TJ but 
for the TV channel. The experimental data are taken from 
Refs. [11, |H. 

RPA based on the covariant density functional theory. 
The fine structure of SD excitations in the most up-to- 
date 16 0(p, n) 16 F experiment is excellently reproduced 
without any readjustment in the functional. The exist- 
ing discrepancy between the data and the shell model 
calculations for the J x = 2 _ "shoulder" structure at 
E x ~ 6 MeV has been clarified. The characteristics of SD 
excitations are understood with the delicate balance be- 
tween the tr- and w-meson fields via the exchange terms 
without the necessity of explicit tensor interactions. The 
fine structure of SD excitations for 16 0(n,p) 16 N channel 
has also been predicted for future experiments. 
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